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Abstract 

We study stochastic properties of the empty space for stationary germ-grain 
models in M. d , in particular we deal with the inner radius of the empty space with 
i -^h ■ respect to a general structuring element which is allowed to be lower-dimensional. 

We consider Poisson cluster germ-grain models and Boolean models with grains that 
are clusters of convex bodies and show that more variable size of clusters results 
in stochastically greater empty space in terms of the empty space hazard function, 
i . We also study impact of clusters being more spread in the space on the value of the 

empty space hazard. Further we obtain asymptotic behavior of the empty space 
' hazard functions at zero and at infinity. 
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: 1 Introduction 

The statistical analysis of a spatial pattern Z C l rf is based on the assumption that Z 
is a random set in M. d . Distance methods for point patterns usually begin by estimating 
nonparametric summary functions such as the empty space function, nearest neighbor dis- 
tance distribution, Ripley's i^-function and derived statistics such as the pair correlation 
and the J-function. For a spatial pattern Z C M. d , which is not necessarily a point process, 
such as germ-grain model a particularly useful functional for estimating properties of Z is 
the empty space hazard (rate) function. Smaller empty space hazard function in a germ- 
grain pattern with the same germ intensity intuitively corresponds to a more clustered 
pattern with more empty space. In this paper we shall give some sufficient conditions 
for the hazard rate ordering of the empty space distribution in two compared germ-grain 
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models using particular classes of such spatial patterns such as Poisson cluster point pro- 
cesses, Poisson cluster germ-grain models, and mixed Poisson germ-grain models. Our 
sufficient conditions will be given in a language of stochastic orderings which we define 
for some parameters of these models. Some early results on stochastic ordering of random 
closed sets can be found in Stoyan and Stoyan [21], and a characterization of the strong 
stochastic ordering of random closed sets is given by Molchanov [15] in Theorem 4.42. 
However, apart from Section 3.8 of Hall [6] (dealing with volume fractions) and Last and 
Holtmann [12] (dealing with the spherical contact distribution of a Gauss- Poisson model) 
we are not aware of papers comparing functionals of stationary random closed sets. We 
introduce two stochastic ordering relations which are useful for comparison two distribu- 
tions with equal expected values in which case both relations imply that a larger variable 
in these orderings has a bigger variance. It turns out that the impact of a larger variance 
in a spatial pattern for some component of a germ grain model with fixed expectation is 
that a larger variance gives more empty space with more clustering. But the situation 
is more complicated if we compare models with equal intensities of Poisson cluster germ 
points but at the same time increasing or decreasing intensities of underlying Poisson 
processes. 

To be more precise, for two random variables rj, fj, taking values 0, 1, 2, . . . we say 
that they are ordered in the length biased probability generating functions ordering, and 
write r) <i~ g fj if the corresponding length biased variables r)i and fji are ordered in 
the probability generating ordering, see [191 Section 1.8]. For two nonnegative random 
variables A 1; A 2 we say that they are ordered in the first cumulant order and write 
A < cum A if for the corresponding cumulant generating functions Ca, C a , taking the first 
derivatives we have C' A (s) > CUs), s G [—1,0]. We show for Neyman-Scott processes with 
cluster sizes r\, fj that r\ <i- g fj implies that the corresponding empty space distributions 
are ordered in the hazard rate ordering. For mixed Poisson germ-grain models we show 
that if the random intensities are ordered A < cum A then the corresponding empty space 
distributions are ordered in the hazard rate ordering. For both introduced orderings, 
under the additional assumption that the ordered variables have equal expected values, 
the greater variable has larger variance. Therefore the mixed Poisson germ-grain model 
behaves similarly as the Neyman-Scott germ-grain model with respect to variance changes 
inside the model. We study also behavior of the Gauss-Poisson model. For details, 
additional properties of the orderings, and examples see Sections 4 and 5. We shall also 
study there the asymptotic behavior of empty space hazard rates at zero and at infinity. 

2 Preliminaries 

2.1 Empty space hazard functions 

Let us recall the definition of the empty space distribution function F (called also the 
first contact distribution function). For ||x|| the Euclidean norm of a vector x G Mr, and 
d(x,A) := inf{||x — y\\ : y G A} the distance between x G M. d and a set A C M. d , this 
distribution is given by 

F(t) = P(d(0,Z) < t), t>0, 
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where denotes the zero vector. The value F(0) = P(0 G Z) is the volume fraction of 
Z. F can also be written in terms of the capacity functional Tz (defined by Tz(K) := 
P(Z n K ^ 0), if C M d compact), namely F(t) = P(Z fl ^ 0) = T z (B(x,t)), 

t > 0, where -B(x, i) is the closed ball with center x and radius t. Stationarity of Z ensures 
that F(t) does not depend on x. Hence, F(t) is the probability that Z hits the ball 
B(t) := B(0,t). There are many reasons for studying other distances than the Euclidean 
distance. For example digital image analysers estimate rather polygonal distance than 
spherical one. To estimate isotropy of point patterns one needs elliptical distance. It is 
clear that the distribution of Z is not determined by Tz(B(x, t)) for all balls, and a larger 
class of sets than the class of balls provides a better information on Z. 

The usual way of introducing other distances than the spherical one is to fix a struc- 
turing element (gauge body) B C M d . This is a compact convex set having G B. Then 
the .B-distance of a point x G lR d to a set K C M d is defined by 

d B (x, K) := inf{* > : (x + tB) fl K ^ 0}. 

It is possible that the set on the right side is empty, e.g. if B is lower-dimensional. In 
such a case we set d B (x, K) := oo. Note that we have the translation invariance property 
d B (z + x, z + K) = d B (x,K), for all z G M. d . Clearly, d B (x,K) < t if and only if x 
is contained in the generalized outer parallel set K + tB* of if, where B* denotes the 
reflected set {— x : x G B}. If B is full dimensional (i.e. has a non-empty interior) and 
centrally symmetric (i.e. B* = B), then d B (-, •) is a metric on IR d induced by the norm 
d B (-,0), and the pair (M d ,(i B ) is called a Minkowski space. 

If the 5-distance ds{x,K) of a point x ^ K is attained in a unique point y in the 
boundary dK of if (that means, if (x + ds(x, K)B) n K — {y}), then we define the 
relative metric projection of x on if by ps{K,x) := y, and the contact direction vector 
u B (K,x) as the element of 95* given by 

tv \ x ~ y 

u B {K,x) := — — — . 

d B (x,K) 

The points x £M. d \K for which the distance d B (x, K) is attained in more than one point 
of K (and for which u B (x,K) is therefore not defined) form the exoskeleton exo B (K) of 
K (see Hug, Last and Weil [IH])- In the Euclidean case, and if if is a finite or locally 
finite set, exo B m(K) is (the boundary of) the Voronoi tessellation generated by K. 
We define the directed, B -relative empty space function F B of Z by 

F B (t, C) := P(d B (0, Z) < t, u B (Z, 0) G C), t > 0, C G B\ (2.1) 

where B d is the system of Borel subsets of M d . Here we use the convention u B (Z, 0) := u 
if G 2 or d B (0,Z) = oo, where Wo £ $-5* is fixed. Definition (12. ip is subject to the 
assumption that the vector u B (Z, 0) is P-a.s. well-defined on {0 < d(Q, Z) < oo}. If Z is a 
random closed set and i? is strictly convex, containing in its interior, then this is indeed 
the case. This follows from the fact that exo B (Z) has volume and from stationarity 
of Z. (More general cases require a suitable assumption on the relative positions of Z 
and B, see Subsection 12.21) . The function F B determines the joint distribution of the 
pair (d B (0, Z),u B (Z, 0)), and hence that of the contact vector d B (0, Z)u B (Z, 0). For each 
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fixed t, the function F B (t,-) is a measure on R d concentrated on dB* . The function 
F B (-) '■= F B (-,M. d ) is called (B-relative) empty space function of Z. 

The Minkowski addition of two sets C, D C M. d (C (B D := {x + y : x <E C,y <E D}) 
gives another form for F B . Stationarity easily implies that 

F B (t) :=F B (t,R d ) =V d (A)- 1 E[V d ((Z®tB*)nA)}, t > 0, (2.2) 

for each Borel test set A, such that the volume Vd(A) of A is positive and finite, see e.g. 

HUE]- 

Hansen, Baddeley and Gill [7], utilizing Federer's coarea theorem, showed that the 
empty space function F B of a random closed set Z is absolutely continuous on (0, oo) 
with density 

f B (t) = V d (AY 1 E \ [ \\Vd B {s,Z)\\- l U d -\ds) 

_JAnd(Z@tB*) 

where % l , % G {0, . . . , d}, denotes i-dimensional Hausdorff measure on M. d , OA denotes the 
boundary of A, and Vd B denotes the gradient of the function d B . In the Euclidean case, 
this formula reduces to 

f{t) = V d (A)- 1 E[H d -\d(Z © B{t)) n A)], (2.4) 

In this case the empty space hazard equals the ratio of the expected measure of the 
boundary d(Z © B(t)) inside the window A to the volume of of the space not occupied 
by Z © B(t) inside A. 

For general B the empty space hazard is given by 



(2.3) 



rsit) V d (A)-E[V d (Z®tB*)nA)] E 



(2.5) 



WVdB&Z^U^ids) 

<And{zmB*) 

which is intuitive in the sense that the empty space hazard depends on the speed of 
increase of the distance function ds along all coordinates. It is possible to rewrite this 
formula in terms of the support function Kb of the gauge body B (see [7j). 



rBit) V d (A) - E[V d (Z © tB*) n A)} E 



[ {h B {uB{Z,s)))- x H d - x {ds) 
JAnd(zmB*) 



■ (2-6) 



U And(Z®tB*) 

The direction dependent (sub) distribution functions, Fb(-,C), as defined in (12. ip are 
also absolutely continuous on (0, oo) for any C G B d . Letting f B (-, C) denote its density, 
we define 

where a/0 := for all a G K. We call the function rs(-,C) the directed, B-relative empty 
space hazard of Z. 

For fixed B, we shall order two random sets Z, and Z with respect to their S-relative 
empty space hazard functions, and write 

Z < h - B Z (2.8) 

iff for all t > 0, and C G B d 

r B {t,C)>f B {t,C), t>0. (2.9) 

This ordering is stronger than the usual (strong) stochastic ordering of empty space 
distributions, for further details on such orderings see e.g. Szekli [22j Section 1.4]. 
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2.2 Empty space hazard rates via support measures 

By a germ- grain model in M. d we mean a random set of the form 

oo oo 

Z= \J(X n + Q= \J{x + £ n :xeX n }, 

n=l n=l 

where the random points £„, iigN, represent the locations of the germs and the primary 
grains X n , n G N, are assumed to be random non-empty compact subsets of M. d . We 
assume that the (simple) point process N : = {£„ : n G N} is stationary, that is the 
distribution of the shifted point process N + x := {£„, + x : n G N} does not depend on 
x G IR d , and that N is independent of (X n ) n >i which is a sequence of independent and 
identically distributed random sets. The intensity A := E ca,id{n G N : („ G [0, l] d } of iV 
is assumed to be finite. An important special case is the Boolean model, where the germs 
are located according to a homogenous Poisson process. The grains X n as well as the 
germ-grain model Z itself are measurable mappings from Q into the set F of all closed 
subsets of M. d . Measurability refers to the smallest a-field of subsets of F, containing the 
sets F^ = {F : F n K ^ 0}, for all compact sets K C M d . Stationarity of N entails 
that also Z is stationary, i.e. that the distribution of Z + x does not depend on x. It is 
convenient to denote by X Q a typical grain having its distribution equal to that of X n . 
We assume that E[Vd(X + K)\ is finite for all compact K C M. d . We will use later that 
the capacity functional of a Boolean model is given by 

P(Z n K ^ 0) = 1 - exp[-XE[V d {X + K*)}}. (2.10) 

In particular, the volume fraction of a Boolean model is given by 

P(0€Z) = 1- exp[-XE[V d (X )}}. (2.11) 

We refer to Stoyan, Kendall and Mecke (20] and Schneider and Weil [H] for a detailed 
introduction to germ-grain models. 

Consider a convex, compact and non-empty set K C M d . We assume that .fT and B* are 
in general relative position, which means that K and B* do not contain parallel segments 
in parallel and equally oriented support (hyper)planes. This means that K and B have 
independent support sets, see [HJ p. 611] for more detail. A sufficient condition is that K 
or B is strictly convex, This assumption guarantees that ps(K,x) (and hence Ub(K, x)) 
is defined for all x ^ K. Then there are finite measures Cq{K; B; •),... , Cd-x{K; B; ■) on 
M. d x M d which satisfy the local Steiner formula 

V d ({x Gl ,i :0< d B (x, K) < t, (p B (K, x),u B {K, x)) G A x C}) 

= t^bd-iQiK; B-AxC) (2.12) 

for all A, (17 G £> d , where 6j (i G N) denotes the volume of the unit ball in W, and 
bo := 1. These relative support measures of K are uniquely determined by (I2.12p . They 
are concentrated on OK x dB* and in fact on the relative normal bundle 

N B (K) := {(p B (K,x),u B (K,x)) : x £ K} 
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of K. If B = -B(l) then the measures Ci(K;-) := Ci(K; B(l); •) are the generalized 
curvature measures of K. The total mass Vi(K) := Ci(K;M. d x M d ) is the ith intrinsic 
volume of K. In particular, V^(i^) is the volume of K, Vd-i(K) equals one half of the 
surface area, Vd-2 (K) is proportional to the integral mean curvature, V\{K) is proportional 
to the mean width of K, and Vq(K) = 1. Equation (12.121) implies the classical Steiner 
formula 

d 

V d (K © B{t)) = b d ^t d -%(K). (2.13) 

In the general case, the total mass Ci(K\ B\ x M d ) is a special mixed volume^ namely 
C i (K;B-R d xR d ) =b^ i r^V(K[i\,B*[d-i}), i = 0,...,d-l. (2.14) 
For % = we have 

C (K; B; R d x R d ) = bfV^B*), (2.15) 

see [18] for the notation used here and for further details on support and curvature mea- 
sures. 

Consider now a germ-grain model Z with convex, compact grains. Assume that the 
reduced second moment measure of iV on IR d , defined by 



E 



i{xe[o,i] d }{x- y e-} 



x,yeN 



(2.16) 



is absolutely continuous, and assume that the typical grain X and B* are a.s. in general 
relative position. It then follows, that Ub(Z, 0) is almost surely well-defined on {0 < 
g?b(0, Z) < oo}. This can be proved as Proposition 4.9 in [9]. 

By the Steiner formula ( 12.131) our general integrability assumption on X (see Subsec- 
tion 12. 2p is equivalent to the finiteness of the mean intrinsic volumes 

Vi:=EVi(X ), i = 0,...,d, (2.17) 

of the typical grain. The Steiner formula ( 12 . 1 3[) together with the local Steiner formula 
f l2~T2|) imply that 

V i>B ■= ECi{X ; B;R d x R d ), i = 0,...,d-l, (2.18) 

are finite as well. Therefore the mean relative support measures of the typical grain, 
defined by 

C iiB (-):=ECi{X ;B;-), % = 0, . . . , d - 1, 

are finite measures on M. d x M. d . 

We further use the Palm probability P 1 ^ of P with respect to N (see (U ED])- We 
can interpret P^(A) as the conditional probability of the event A e T given that is 
a "randomly chosen point" of N. Let us define X(x) := X n if x = ^ n for some n, and 
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X(x) := 0, otherwise. Then under P^, {(x,X(x)) : x G N} is an independently marked 
point process, i.e. conditionally on N, the grains {X(x) : x G N} are independent and 
have the same distribution as Xq. For V^b > 0, t > 0, and C G £> d , let 

G ijB (t, C) := J l{u G C}P° N (d B (y + tu, Z l ) < t)C iiB (d(y, u)), (2.19) 

where 

Z--.= (J (X(x)+x) (2.20) 

rreJV\{0} 

is the union of all grains except for the grain located at the origin. We set G^b = for 
Vi t B = 0. The function Gj ; B(-,lR d ) can be interpreted as the distribution function of a 
random variable £, say, which can be constructed as follows. First one selects a point Y of 
N at random. Then one samples a random element (X, W) according to the distribution 
V^C^b- If Y is not covered by U xe N\{Y}(X(x) +x), then £ is the £>-distance from Y + X 
to the exoskeleton exos{Z), in the direction W. Otherwise £ = 0. We shall utilize the 
following functions 

JMt.c)-.= GiAc °: c) :% B(t - c) , i = o,...,d-i, (2.2i) 

where z G {0, . . . , d - 1}, t > 0, C C R d x M d is a Borel set, and 

C iiB (oo,C) := lim G hB (t,C) = Vr'C hB (R d x C). 

t—too ' 

Special cases of these functions were introduced in (T2J Section 5] after the point process 
case had been treated in [Ti] . The functions Ji, B (t, C) can be used as non-parametric mea- 
sures for expressing differences between a general germ-grain model and Boolean model 
with the same values of XVi jB . Intuitively speaking, such measures detect interactions 
and clustering effects. In the Euclidean case (and for C = M. d ) the following theorem was 
proved in [HI [12]. The following relative version is implicit in Hug and Last [9], at least 
in the case of a strictly convex B. The general result can be derived from Theorem 5.1 
in Hug, Last and Weil pEj. 

Theorem 2.1. Consider a stationary germ- grain model satisfying the assumptions for- 
mulated above. Then F B is absolutely continuous and the B-relative empty space hazard 
r B is given by 

r B (t,C) = Ysid-i^-X-AV^J^C). (2.22) 

If N is a Poisson process (i.e. Z is the Boolean model with convex grains), then 
Slivnyak's theorem (see e.g. Stoyan, Kendall and Mecke [20]) implies that 

Pn(Mv + tu > Z-)>t) = l- F B {t). (2.23) 

Hence f)2.22p simplifies to 

d-l 

r B (t, C) = J^id - t^t^XSi^C), (2.24) 

i=0 
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where S^b(C) := C , j ) s(]R d x C). In the case of a strictly convex gauge body B this 
result can be found in [9]. Note that in the Boolean model rg(-, M, d ) is determined by the 
intensity A and the mean mixed volumes V^b, ■ ■ ■ , V d -i,B of X . For d = 2 and 5 = B(l), 
for instance, the only parameter of Xq influencing the empty space hazard rate is its mean 
boundary length V±. Note that in the Boolean model with convex grains the empty space 
hazard rate is increasing, and asymptotically 

lim t^rsit, C) = \db d S , B {C) = Xdb d E[C {X ; B; R d x C)]. 

t—^oo 

By ( I2.15P we obtain in particular that 

lim * 1 - < V B (t > R << ) = \db d E[b^ x V d {B*)\ = XdV d {B*). 

t— Kxi 



3 Results for Poisson cluster point processes 

In this section we shall consider germ-grain models where the germ process will be a 
Poisson cluster point process, and grains will be one point grains attached to germs, 
hence Z = N. Such a model is a special case of a germ-grain model with convex (one 
point) grains but, alternatively, this model might be seen as a Boolean model (germs form 
a Poisson process) with non-convex grains (point clusters). Recall from j3] that a Poisson 
cluster point process can be written as 

N={jL x + x, (3.1) 

where II is a Poisson process with positive and finite intensity An and the family {L x : 
x G 11} consists of finite random point processes on M. d . Given II, the family {L x : x G 11} 
is i.i.d. with the same distribution as a typical cluster L . We assume that 

7 := E card L 

is finite and positive, hence A = An7- 
Example 3.1. Assume that 

io=| 9 ' , if " = °' (3.2) 

[{y ijn :i = l,...,n}, if»7 = 7i>l, 

where the random cluster size rj > 0, and the random vectors Yj >n , n G N, i = 1, . . . , n, 
are independent. Assume also that the Yj jn have the same (cluster point) distribution V, 
say. Then N is called Neyman- Scott process. We always assume that 7 = E\rf[ is positive 
and finite. 

Example 3.2. Let rj be a {l,2}-valued random variable and assume that 

{0} ' ifr? = 1 ' (3 3) 

{0,F}, if rj = 2, 

where Y is a random vector independent of 77. Then N is called Gauss-Poisson process. 
In this case a cluster L + x associated with a parent point x say, contains x and, with 
probability p := P(r] = 2), also a secondary point x. Note that the mean cluster size is 
given by 7 = 1 + p. 
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In this paper we always assume that the reduced second moment measure 



E 



l{x-ye-} 



x,y£L 



(3.4) 



of L is absolutely continuous on ~R d . The well-known second order properties of Pois- 
son cluster point processes (see e.g. [I]) easily imply that the measure defined at (12.161) 
is absolutely continuous as well. For the Neyman-Scott process of Example 13.11 our as- 
sumption on L is implied by the absolute continuity of the cluster point distribution V. 
For a Gauss-Poisson process it is sufficient to assume that the seondary point Y has an 
absolutely continuous distribution. 

In this section we fix a gauge body B that contains a non-empty neighborhood of 0. 
This way we exclude the trivial case F B (t) — 0, t > 0. Let v B be the measure on IR d given 
by 



d / l{x/d B (Q, x) £ -}dx, 



(3.5) 



where d B (0,x) := 6^(0, {x}). Using this measure we can express the empty space hazard 
of TV as follows. 

Proposition 3.3. The B-relative empty space hazard of the Poisson cluster point process 
Z = N is given by 



v B (du). (3.6) 



r B (t, C) = An^" 1 I E f 1{((L - x) \ {0}) n (tu + tB) = ®}L {dx) 

JC ]_J9. d 

Proof. Computing the left-hand side of (I2.12p for K = {0} and A = {0} x C easily 
shows that 

db d C ({0};B;{0}x-) = u B (-) 

and that Cj({0}; B; •) = for i > 1. The result is then a consequence of Proposition 14.11 
below. □ 

If Lq = {0} we have N = U and 

r B (t,C) = Xt d - 1 u B (C), (3.7) 

in accordance with (12.241) . 

If B = B(l) then u B is the (d — l)-dimensional Hausdorff measure on the unit sphere 
S d-i ._ g B ^ and 



r B (i)(t,C) = A n t 



d-l 



E 



c 



1{((L - x) \ {0}) n B(tu, t) = ®}L (dx) 



H d -\du). 

(3.8) 

Example 3.4. Assume that iV is a Neyman-Scott process as defined in Example 13. 1[ 
From (13. 6p and a straightforward calculation, 



r B (t, C) = Ant d-1 / / g'(P{Y 1A - x £ tu + tB))V {dx)v B {du) 



(3.9) 



where g' is the derivative of the probability generating function g of rj. This result gener- 
alizes formula (30) in [TT] . 
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Example 3.5. Assume that N is a Gauss-Poisson process as defined in Example 13.51 
Then 



(An) 



-i t i- d 



r B {t,C) ={l-p)v B {C) (3.10) 

+ P P(Y <£tu + tB)v B (du)+p / P(-Y <£tu + tB)v B (du). 
Jc Jc 



Let N be a general Poisson cluster point process as above. It is interesting to note 
that t 1 ~ d r B (t, C) is monotone decreasing in t. This is a direct consequence of (13. 6p . Then 
it is instructive to compare r B with the right side of (13 ,7p . From monotone convergence 



limt 

t-»o 



l-d 



r B (t,C) = \du B (C), 



(3.11) 



i.e. for small values of t the empty space hazard of a Poisson cluster point process behaves 
approximately like the empty space hazard of a Poisson process with the same intensity. 

Next we deal with the asymptotics of the empty space hazard for large t. It turns out 
that it is the same as that of the Poisson process II thinned at points x where L x is empty, 
for which intensity equals P(L ^ 0)An- This means in a sense, that points in clusters 
cannot be distinguished from a very far distance, irrespective of any specific assumptions 
on the typical cluster Lq. This is generalizing equation (22) in [TT| . A weaker version of 
this latter result has been rediscovered by Bordenave and Torrisi [3]. 

Proposition 3.6. The B-relative empty space hazard of a Poisson cluster point process 
satisfies 



lim t l - d r B {t,C) = P(L jt 0)A n ^(C). 

t— >oo 



(3.12) 



Proof. The tangential cone (or support cone) T(B,u) of B at u G B is the closure 
of T'(B,u) := {t(x — u) : t > 0,x E B}, see Schneider [T7j. From (I3.6P and monotone 



convergence 



lim t l ~ d r B {t, C) = An / E [ 1{((L - x) \ {0}) n T'{B, -u) = 0}L o (c/x) 

Jc U 



u B (du). 
(3.13) 



Since the measure (13 .4p is absolutely continuous it follows that Yq a.s. does not intersect 
the boundary of T(B, —u). Hence 



lim t l - d r B {t, C) = An / E [ 1{((L - x) \ {0}) n T{B, -u) = <&}L (dx) 
*- ) ' 00 Jc U 



v B {du). 
(3.14) 



Let us now fix for a moment a regular boundary point u of B* . This means that B has a 
unique supporting hyperplane at — u, see Schneider [TTJ . Then T(B, —u) is the supporting 
half-space of B at — u, see Section 2.2 in Schneider [17]. It now follows as in [111 Section 
6.5] that 



E 



J 1{((L - x) \ {0}) H T(B, -u) = ®}L (dx) 



P(W0)- 
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It remains to note that z/g-a.a. u £ M d are regular boundary points of B* . This follows 
from the fact that Vb is absolutely continuous with respect to (d— l)-dimensional Hausdorff 
measure on dB* and Theorem 2.2.4 in Schneider [17]. □ 

We shall now consider two Poisson cluster processes N and N having the same inten- 
sity, satisfying the technical assumption formulated at ( 13. 4p and with 5-relative hazard 
rates r B and r B , respectively. We will establish several sets of assumptions implying the 
hazard rate ordering ( 12. 8p . In all cases N will have more clustered patterns in a sense. 
Therefore one may expect more empty space in the Poisson cluster process N. Still it is 
somewhat surprising that this happens in the strong sense of ( I3.17p . In Section H] we will 
state the corresponding results for Poisson cluster grain models. Our first ordering result, 
dealing with Neyman Scott processes, is a special case of Proposition 14.91 which will be 
proved later. 

Let 7/ and fj be two counting variables (i.e. taking values in {0, 1,2,.. .}) and let r\i be 
the (shifted) length-biasedversion of rj. This means that rji has distribution Elrj]^ 1 E[i]l{r]— 
1 £ ■}]. Denoting by fji the length-biased version of fj, we define the length biased proba- 
bility generating functions ordering 

V<i-gV (3-15) 

by 

E[s m ] >E[s% s£[0,l]. 

This means that r)i is smaller than fji in the the generating function order (see Stoyan 
|19j . Section 1.8.) Note that rj <i~ g fj is equivalent to 

E[rj\~ l E[l{q > ljr/a^ 1 ] > E[fj)- 1 E[l{fj > a £ [0, 1]. 

Another way of expressing this relation is by the generating functions of r\ and fj. Denote 
by 

g v {s) :=E[s% sE [0,1], 
the probability generating function of rj. Then (13. 15ft means that 

Ein)- 1 ^) > E^g'-is), s £ [0, 1]. (3.16) 

The relation (13. 15ft does not imply that the corresponding expected values i?^],^^] 
are ordered. For example if 77 = 1, 771 = 2, and 772 equals with probability 1/2 and 
1 with probability 1/2 then 77 <i~ g r/i, and r\ <i- g but E[rf\ = 1 < E[r]i\ = 2, and 
E[r]} = 1 > E[t]2] = 1/2. However, if E\rf[ = E[fj] then (I3.15P implies that Var[?7] < Var[r/], 
therefore this relation is a variability ordering. 

Proposition 3.7. Assume that N and N are Neyman- Scott processes with cluster sizes 
77 and fj, respectively, having the same distribution V of cluster points. If N and N have 
the same intensity and 

V <l-g V 

then 

Z < h - B Z (3.17) 
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Example 3.8. Let iV and N be as in Proposition 13.71 and assume that rj = 1. That is N 
is again a Poisson process with intensity 1 since the points of the original Poisson process 
IT are independently shifted. Let a G [0, 1]. The length biased variable rji = 0, and its 
generating function S[s' B ] = 1, therefore for each fj we have 1 = rj <i~ g fj, and Z <h-B Z. 
This example shows that within the class of Neyman-Scott processes with fixed intensity 
the stochastically smallest empty space appears for pure Poisson germ processes. 

Example 3.9. Let N and iV be as in Proposition 13.71 and assume that rj and fj are 
Poisson distributed with parameter c and c, respectively. (Then the clusters are finite 
Poisson processes.) We have for any a G [0, 1] that 

E[ti\- x E[1{ti > l}^- 1 ] = e~ c(1 - a) . 

Therefore (13. 15[) (and hence Z <h-B Z) holds iff c < c. Note that for Poisson distributed 
rj the length-biased rji have the same (Poisson) distribution. Proportionally more Poisson 
cluster points in a Neyman-Scott process indeed generate more clustering and lead to a 
stochastically larger empty space. 

Example 3.10. Assume that rj and fj are binomial distributed with parameter (n,p) and 
(n,p), respectively. Since g n (s) = ((1 — p) +ps) n , it follows that (13.151) is equivalent to 

((1 -p)+ ps) n - 1 > ((1 -p)+ psf-\ s G [0, 1]. (3.18) 

If, for instance n = n, then this inequality is implied by p < p. 

Example 3.11. Assume that r\ and fj are negative binomial distributed with parameter 
(p,r) and (p,f), respectively. The corresponding length-biased variables rji , fji are again 
negative binomial distributed with parameters (p,r + 1) and (p,f + 1), respectively. For 
p = p, if r < f then 77 <i~ g fj, and for r = f, if p > p then 77 <i- g fj- This is a special case 
of a more general setting. If 77 = Yli=i an d V — Y^i=i ®ii ^ OY nc ^ variables {i?j}i>i, and 
independent of k (all variables taking on natural values) then f)i <i- g $i, and k <i- g R 
implies r\ <i- g fj. 

The next result is a special case of Proposition 14.101 

Proposition 3.12. Assume that N and N are Gauss-Poisson processes having the same 
distribution V of the secondary point. Assume that N and N have the same intensity and 
that the probabilities p and p for having a secondary point satisfy p < p. Then Z <h-b Z . 

In our next result we will multiply each point of the typical cluster Lq of with a 
random variable W < 1, i.e. we assume that the typical cluster L of has the distribution 
of WLq = {Wx : x G L }. Compared with L , the points of L are more spread out. 
Note that we allow any sort of dependence between L and D. 

Proposition 3.13. Consider two Poisson cluster processes N and N based on the same 
Poisson process U and typical clusters L and L , respectively. Assume that L is dis- 
tributed as WLq for some random variable W < 1. Then Z <h-B Z. 
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Proof. Take x±, . . . , x n G M d and set ip := {xj : j — 1, . . . , n}. Let ?i> < 1 and define 
if) := {wa;-,- : j = 1, . . . , n}. Let u G 95* and t > 0. In view of (13.61) it is sufficient to show 
that 

l{((if) - kji.) \ {0}) n (tu + tB) = 0} < 1{((V> - Xi) \ {0}) n (tu + tB) = 0} (3.19) 

holds for any % G {1, . . . ,n}. Assume that the left-hand side of (13.191) equals 1. This is 
equivalent to w(xj — xi) tu+tB for all j ^ i. Since u G dB* we have that G d(tu + tB). 
Therefore the convexity of tu + tB and w < 1 imply that also Xj — Xi ^ tu + tB for all 
j 7^ i. This shows (I3.19P and hence the proposition. □ 

Recall that a random variable W is stochastically smaller than another r.v. W, if 
> t) < P{W > t) for all t > 0. 

Proposition 3.14. Let N and N be as in Proposition \3.1HA Let L be a point process 
and assume that L is distributed as WL for a random variable W > independent of L. 
Assume that L is distributed as WL for a random variable W > independent of L. If 
W is stochastically smaller than W then Z <h-B Z . 

Proof. By inverse coupling based on an uniformly distributed random variable that 
it independent of L we can assume that (W, W) is independent of L and that W < W 
everywhere on the underlying probability space. Since Lq = W/WLq, we can apply 
Proposition 13. 131 with W replaced by W/W < 1. □ 



4 Results for Poisson cluster germ-grain models 



We now return to a general case of germ-grain models where N is a Poisson cluster point 
processes as in the previous section, and grains X n , located at all points of N, are compact 
and convex. In the next two sections we fix a structuring element B such that our general 
assumption made after (12 . 1 6[) holds, that is X\ and B* are a.s. in general relative position. 
Alternatively, we shall treat such a process as a Boolean model with non-convex grains 
which are 

\J(X(y)+y), xeU. 

y&L x 

For a given finite point pattern if) = {x n : n = l,...,m}, denote by T(ip, •) the 
distribution of the random closed set U™ =1 (A n + x n ), where X\, . . . , X m are independent 
with distribution of Xq. We also set T(0, •) := 5%. Let Yq be a random closed set with 
distribution 



m(0 



-E 
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j HAe-}T((L o -y)\{0},dA) 



(4.1) 



This probability measure is describing the distribution of the germ-grain model associ- 
ated with a typical cluster as seen from a randomly chosen cluster point, after removing 
the grain around the chosen point. We assume that Yq and the typical grain X are 
independent and define for i G {0, . . . , d — 1}, t > 0, and a Borel set C <zM. d 



Ki, B (t,C) :-- 



Vi. 



-E 



l{d B (x + tu, Y ) > t}l{u G C}Ci(X ; B; d(x, u)) 



(4.2) 
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The following proposition yields again an explicit formula for the empty space hazard of 
a Poisson cluster germ-grain model, this time in terms of X and Y describing locally a 
cluster. 

Proposition 4.1. The B-relative empty space hazard of a Poisson cluster germ-grain 
model with compact, convex grains is given by 



i-i 



r B (t,C) = ^(d-i^-X-iXV^K^C). 



(4.3) 



i=0 



Proof. Our aim is to use (12.221) . To do so we recall that the Palm probability measure 
of a Poisson cluster process N satisfies 



where 



P° N (N e-) = E 



QU)--=r l E 



l{iVU^G-}Q° (#) 



(4.4) 



4/eLo 

is the Palm distribution of the typical cluster Lo, see e.g. Stoyan, Kendall and Mecke [20] 
As in Last and Holtmann [12], this implies that 

P° N (d B (x + tu, Z-) >*) = (!- F B {t)) J I \{d B {x + tu, A) > t}Y{^ \ {0}, dA)Q° Lo (dijj) . 



By definition (12.191) and definition of Y this means that 
G iiB (oc,C)-G iiB (t,C) 1 



l-F B (t) 



V, 



-E 



i.B 



l{d B (x + tu, Y ) > t]l{u G C}C itB (d(x, u)) 



By definition (14. 2p . the above right-hand side equals K iyB (t,C). Inserting this into (I2.22p 
gives the asserted equation (14. 3p . □ 

Our next proposition deals with the asymptotic behavior of r B (t,C) as t — > or 
t — > oo, respectively. Note that t l ~ d r B (t,C) is monotone decreasing in t. This is a direct 
consequence of (14. 2p . Recall the definition of the tangential cone given in the proof of 
Proposition 13.61 

Proposition 4.2. The B-relative empty space hazard r B (t, C) of a Poisson cluster germ- 
grain model with compact, convex grains satisfies 



\imr B (t,C) = 2X U E 

t— 5-0 



lim t^TB^C) = db d X u E 

t— >oo 



l{x^Y }C d ^(X ;B;d(x,u)) 



>xC 



(4.5) 



/ l{(x + T(B,-u))nY = i&}C (X ;B;d(x,u)) . 

JR d xC 



Proof. The proof is similar to the proof of (13. lip and (I3.14p . 



(4.6) 
□ 
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Remark 4.3. One might wonder whether the right-hand side of (14. 6p is positive for 
C = ~R d . A simple sufficient condition is that P(cardL = 1) > 0, because in this case 
Y"o is with positive probability empty. Another sufficient condition is to assume that B is 
smooth (any boundary point has a unique supporting hyperplane) and that the diameter 
of the typical grain can take arbitrary small positive values with positive probability. 
These assumptions would allow to apply the method of [HI Section 6.5] on a set of 
positive probability. We do not go into further details. 

Remark 4.4. The second assertion of the previous proposition shows in particular that 
F B is light-tailed, i.e. has a finite exponential moment. 

Now we shall compare the P-relative empty space hazard r B of Z with that of another 
Poisson cluster germ-grain model Z with the same intensity A of germs and the same 
typical grain X . Both underlying germ processes are assumed to satisfy the technical 
assumption formulated at (13 .4p . We denote the characteristics of Z by IT, L , K i B , r B 
etc. We begin with a direct consequence of Theorem 14.11 

Proposition 4.5. Let C G B d and t > such that 

K itB (t,C)>k itB (t,C). (4.7) 

Then (fjOj]) holds. 

An immediate consequence of the above proposition is that the relative empty space 
hazard of a Boolean model is always greater than that of a Poisson cluster germ-grain 
model having the same germ intensity and the same typical grain: 

Corollary 4.6. Assume that Z is a Boolean model with typical convex, compact grains 
distributed as X Q , and Z a Poisson cluster germ-grain model with equal intensity, and 
also typical grains distributed as Xq. Then Z <h-B Z . 

Proof. Since Z is a Boolean model, we have Lq = {0} and 7 = 1. Hence, if V^ B > 0, 
Ki, B (t, C) = Vr^C ijB (R d x C) > K itB (t, C). □ 

While a Boolean model has stochastically smaller empty space than a related Poisson 
cluster germ-grain model, it has a greater volume fraction. Under a different set of 
assumptions (more specific Poisson-cluster processes and deterministic but possibly non- 
convex grains) the result was proved in Section 3.8 of [5J. 

Proposition 4.7. Under the assumptions of Corollary \4- 6] 

P(0 G Z) > P(0 G Z). (4.8) 

Proof. The volume fraction of Z is given by (12. lip . On the other hand, Z is also a 
Boolean model, but based on the Poisson process II and with typical (possibly non-convex) 
grain 

X := |J X(x) + x, 



15 



where L is the typical cluster associated with Z and, given L , the family {X(x) : x G Lq} 
consists of independent random closed sets with the same distribution as X . Therefore 
we obtain from (12. lip 

P(0eZ) = l- exp[-X n EV d (X )]. (4.9) 

We have 

EV d (X Q ) <EJ2 Vd(X(x)) = EJ2 E[V d (X(x))\L ] = >yV d . 

Inserting this into (14.91) and comparing with (12. lip , yields the assertion. □ 

Remark 4.8. Consider the hypothesis of Proposition 14.71 The proof of this proposition 
shows that we have equality in (14.81) iff 

V d (X(x)nX(y)) = x,yeL ,x^y P - a.s. (4.10) 

This is, for instance, the case if the cluster points have minimal distance 2t from each 
other for some to > 0, and X is a.s. contained in the ball B(to). 

Next we prove a more general germ-grain version of Proposition 13.71 

Proposition 4.9. Consider two germ-grain models with the same typical grains Xq such 
that N and N are Neyman-Scott processes with cluster sizes i] and fj, respectively, and 
the same cluster point distribution V. If N and N have the same intensity and 

V <i-g V (4.11) 

then Z <h-B Z . 

Proof. By Proposition 14.51 it suffices to show 

P(d B (x + tu, Y ) >t)> P(d B (x + tu, Y ) > t) (4.12) 

for all i,m6K <! and all t > 0. Letting B' := tB + x + tu and using the definition (14. ip of 
the distribution of Y , we obtain that 

E[ V }P(Y n B' = 0) = E ii A n B ' = > r (( L o - v) \ W' dA ) 

oo „ 

= P( V = 1) + P (V = n)n / fivT^Vidy), 

n=2 

where f(y) := J P((Xq + z — y) DB' — $)V(dz) and where the second identity comes from 
a straightforward calculation using the definition of the typical cluster L of a Neyman- 
Scott process. By Fubini's theorem this means that 

E[r,]P{Y n B' = 0) = J Efrfiy^Vidy). 

We can now use our assumption (I3.15P to derive 

P(Y n B' = 0) > Efr]- 1 J E^fiyf-^Vidy). 

Reversing the above steps, we get (14.121) and hence the asserted result. □ 
Our next result generalizes Theorem 5.4 in |12j . 
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Proposition 4.10. Consider two germ-grain models with the same typical grains X based 
on Gauss-Poisson processes N and N . Assume that N and N have the same intensity and 
that the probabilities p and p for having a secondary point satisfy p < p. Then Z <h-B Z . 

Proof. Fix t > and C G B d . From the defining properties of a Gauss-Poisson process 
(see Example I3.5f) we have 

V iiB K ijB (t, C) =\^C hB (R d xC) + [ a(x, u)C hB (d(x, it)) 

p f 

+ — — / b(x,u)Ci )B (d(x,u)), 

- 1 - + P JR d xC 

where 

a(x, u) := P(d B (x + tu,Y + X ) > t), b(x, u) := P(d B (x + tu, -Y + X ) > t), 

where Y and Xo are independent. Therefore, by Proposition 14.51 it suffices to show that 

I — p pa pb 1 — p pa pb 
1 + p 1 + p 1 + p ~ 1 + p 1 +p 1 +p' 

for all a, b G [0, 1]. Simple algebra shows that this inequality is equivalent to 

2p — ap — bp > 2p — ap — bp. 
The latter is implied by our assumption p < p. □ 

5 Results for Mixed Poisson germ-grain models 

In this section we consider a germ-grain model Z based on a mixed Poisson process N. 
This means that there is a random variable A > such that the conditional distribution 
of N given A is that of a stationary Poisson process with intensity A. We assume that 
E[A] (the intensity of N) is positive and finite. 
It is convenient to use the notation 

H B (t) = E\V d (X +tB% t>0. 

Proposition 5.1. The B-relative empty space hazard of a mixed Poisson germ-grain 
model with compact, convex grains is given by 

d-l 

r B (t,C) = ^(rf- z)t d - J - 1 6^E[exp[-AiJ B (t)]]- 1 J E[Aexp[-Aif B (t)]]^ B (t,C). (5.1) 

i=0 

Proof. Again we will use (12.221) . To do so, we note that the Palm probability measure 
PJv of a mixed Poisson process N satisfies 

P° N ((A, AO G •) = £ [Al-^AlftA, N U {0}) G ■}]. (5.2) 

This formula can be derived by conditioning and using the properties of a Poisson process. 
Since, moreover, the conditional distribution P%(Z [ G -|A) (cf. (12.201) for the definition of 
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the random set Z l ) is that of a Boolean model with germ intensity A, we obtain for all 

(x, u) e R d x R d that 

Pn(cIb(x + tu, Z-) >t) = E[A\- 1 E[KP(d B {x + tu,Z)>t\ A)] 

= PA^PAP^O, Z)>t\ A)] 
= P[A]- 1 P[Aexp[-A^(t)]], 

where we have used (12.101) to obtain the last identity. Again by conditioning and (I2.10p we 
have that 1 — Pb(£) = P[exp[— Aife(t)]]. Inserting our findings into the general formula 
(1Z22|) yields the assertion (jOl . □ 

In order to state some stochastic ordering consequences of Proposition EH we introduce 
a stochastic order using cumulants. For two nonnegative random variables A, A we 
say that they are ordered in the first cumulant order and write A < cum A if for the 
corresponding cumulant generating functions C\, C A , taking the first derivatives we have 
C' A (s) > CUs), s G [—1,0]. Note that A < cum A is equivalent to 

£?[exp[-Aa]] _1 E[Aexp[-As]] > Efexpf-As^-^Aexpf-As]], s > 0. (5.3) 

The left-hand side of ( 15. 3 j) is the logarithmic derivative of the Laplace transform s h- > 
P[exp[— As]]. It is also the hazard rate of the distribution function G\, defined by 

G A (s) := 1 - E[exp[-As]], s > 0. 

This is a mixture of exponential distributions. Equation ( 15. 3 j) then means that the cor- 
responding variables are ordered in the hazard rate order, i.e. Ga <h G A . Note that 
for A, A with equal expected values A < cum A implies that Var[A] < Var[A], therefore, 
similarly to the relation <i- g , the relation < cum is a variability ordering in the case of a 
fixed expectation. 

An immediate consequence of Proposition 15.11 is the following counterpart of Proposi- 
tions 14.91 and 14.101 We use similar notation. Again, intuitively speaking, more variability 
in the mixed Poisson model results in a stochastically greater empty space (a stochastically 
larger clustering). 

Proposition 5.2. Consider two germ-grain models with the same typical grains Xq based 
on mixed Poisson processes N and N with random intensities A and A, respectively. 
Assume that A < cum A then Z <h-B Z . 

Example 5.3. Let and A^ be as in Proposition 15.11 and assume that A is Gamma 
distributed with shape and scale parameter a > and /3 > 0, respectively. This means 
that A has density (3°T(a)~ 1 x a ~ 1 exp[—/3x]. The Laplace transform of A can be computed 

as 

ga 

E[exp[-sA}} 



while an equally easy calculation gives 

a (3 a 



£[Aexp[-sA]] 



f3 + s(f3 + sY 
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Assume now that A is Gamma distributed with parameters a and {3, respectively. Then 
assumption (I5.3P means that a/(/3 + s) > a/($ + s) holds for all s > 0. This is equivalent 
to 

*>a, ->,. (5.4) 

Depending on whether or not /? > /3, only one of these equations is relevant. By Proposi- 
tion (15. 4p implies the empty space hazard ordering (13.171) . Assume for instance that 
A is exponentially distributed with mean 1, i.e. a — (3 — 1 and assume furthermore that 
A has also mean 1, that is a — /3. Then (15 .4p is equivalent to a > a. Note that the 
variance of A satisfies Var[A] = a//3 2 = 1//3 < Var[A] = 1 if a > 1. 

As in the Poisson cluster case it follows that the relative empty space hazard of a 
Boolean model is greater than that of a mixed Poisson germ-grain model with the same 
germ intensity. 

Corollary 5.4. Assume that Z is a Boolean model with typical convex, compact grains 
distributed as X Q , and Z a mixed Poisson germ-grain model with equal intensity, and also 
typical grains distributed as X . Then Z <u-b Z. 

Proof. Let A denote the germ intensity of the Boolean model Z and let A be the random 
intensity of the mixed Poisson process underlying Z. It is assumed that E[A] = A. We 
check that condition (15. 3p holds with A = A. This condition means that 

A£[exp[-As]] > £[Aexp[-As]], s > 0. 

In other words: the covariance between A and — exp[— As] has to be non-negative. This 
fact follows from a very well-known statement that a single random variable is associated, 
see Esary et al. [5]. □ 

For completeness we provide the mixed Poisson analogue of Corollary 14.71 The result 
can be found in Section 3.8 of [6] for the more general case of stationary Cox processes with 
an absolutely continuos intensity measure. Although our proof below can be extended to 
arbitrary stationary Cox processes we stick to the mixed Poisson case. 



Proposition 5.5. Under the assumptions of Corollary \5. 4 



P(0 G Z) > P(0 G Z). (5.5) 

Proof. By conditioning and (12. lip . 

l-P(0eZ) = E[exp[-AE[V d (X )}}}. 

By Jensen's inequality this is bounded from below by exp[— A-EfV^Ao)]]. This lower 
bound is just 1 — P(0 G Z). □ 
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6 Concluding remarks 



We have derived several variability properties of the empty space function of Poisson 
cluster and mixed Poisson germ-grain models. It would be worthwhile to study also other 
classes of germ processes. Another interesting task is to find a good notion of spread out 
for a finite point process (with respect to the origin). Proposition 13.121 and Proposition 
!3.13l should be both special cases of the same principle. The first proposition is generalized 
by Proposition 14.101 We believe that also Proposition 13. 131 has a germ-grain counterpart. 

In this paper we have always fixed the distribution of the typical grain. However, it 
would be quite interesting to study the variability of empty space in germ-grain models 
for a fixed germ-process but variable grain distribution. For instance, one might compare 
models with equal expected volumes of the typical grains. To illustrate this task we give 
one simple example that is closely related to some of the results in [21] . 

Example 6.1. Let X be a random convex body such that E[Vd(Xo + K)] is finite for all 
compact K C M d . Let R and R be positive random variables with a finite d-th moment 
and assume that X and R (resp. X and R) are independent. Consider two Boolean 
model Z and Z based on the same Poisson process N and typical grains RX and RX , 
respectively. If 

E[R l ]>E[R% i = l,...,d-l, 

then (12.81) holds for all structuring elements B such that Xo and B* are a.s. in general 
relative position, and all Borel sets C C M. d . This follows from ( 12. 24ft and the scaling 
property 

Ci(aX ; B; R d x C) = dC^X^ B; R d x C), a > 0, 
see e.g. [T7] for the Euclidean case B = 5(1). 
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